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ABSTRACT 

Among the biggest surprises revealed by COBE and confirmed by WMAP measurements of the 
temperature anisotropy of the CMB are the anomalous features in the 2-point angular correlation 
function on very large angular scales. In particular, the i = 2 quadrupole and t = 3 octopolc 
terms are surprisingly planar and aligned with one another, which is highly unlikely for a statistically 
isotropic Gaussian random field, and the axis of the combined low-<? signal is perpendicular to the 
ecliptic plane and the plane defined by the dipole direction. Although this < 0.1% 3-axis alignment 
might be explained as a statistical fluke, it is certainly an uncomfortable one, which has prompted 
numerous exotic explanations as well as the now well known "Axis of Evil" (AOE) nickname. Here, 
we present a novel explanation for the AOE as the result of weak lensing of the CMB dipole by large 
scale structures in the local universe, and demonstrate that the effect is qualitatively correct and of a 
magnitude sufficient to fully explain the anomaly. 

Subject headings: cosmology: Lensing — cosmology: large-scale structure 



1. INTRODUCTION 

The full-sky maps obtained by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) in its first year 
of observ ation have revolutionized the study of the 
CMB skv HBermett et «,1 jlSfitfiafl IHinshaw et a,l]l2003at 
ISpergel et al.l l2003). One of the many exciting uses of 
these maps is the window they offer u s on the physics of 
the early universe ijPeiris et al.l l2003'l and the standard 
inflationary model. This picture predicts a statistically 
isotropic Gaussian random CMB temperature anisotropy 
to an excellent approximation, so that WMAP has en- 
abled considerable study of non-Gauss ianity and statis- 
tical anisotropy of the CMB (see e.g. ICopi et alJ 120051 
and references therein). 

Although WMAP data has largely been consistent with 
the standard cosmological picture, some challenges to the 
standard model have emerged. Among the most interest- 
ing occur on the largest angular scales, and arc therefore 
called the "low-£ anomalies" . The first of these, the sur- 
prisingly small quadru pole moment, was not ed a decade 
ag( > in the COBE data ( Hinsha w et al.lll996|) and has re - 
cently been confirmed bv WMAP ijSpergel et al.ll2003|) . 
In addition, it is now known that the oct opole is highly 
planar and aligned with the quadr upole ijSchwarz et alJ 
120041: Ide Oliveira-Costa et al.11200^ . and that the planes 
defined by these multi-poles is perpendicular to both the 
dipole and the ecliptic. 

The original low quadrupole anomaly has long been 
dismissed as either the result of some residual systematic 
error or as a statistical fluke. However, the higher quality 
of data now available from WMA P strongly challenge th e 
residual systematic explanation ijTegmark et al] [20031. 
and while a statistical fluke cannot be ruled out, the odds 
against are unco mfortably long. In one recent study, 
ICopi et all lj2005T) have used the multi-pole vector formal- 
ism to show that a purely accidental alignment is unlikely 
in excess of 99.9%. They have also shown that most of 
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the 1 = 1 and i = 3 multi-pole vectors of known Galactic 
foregrounds are located far away from those observed in 
WMAP data, strongly suggesting that residual contam- 
ination by foregrounds which are currently included in 
the analysis is not a viable explanation. It is precisely 
this combination of a complete lack of any known sys- 
tematic error, and long odds against random alignment 
that has earned the low-£ ali gnment anomaly the nick - 
name "Axis of Evil" fAOE: lLand fc Magueiid I2005a|) . 
The axis appears to possess mirror symmetric proper- 
ties ([Land fe Magueiial2005bj) and may be related to the 
north-south asymm etry in the angular power spectrum 
(Erik sen et alll200l . 

The fact that the AOE alignment includes sources of 
both cosmological and local origin suggests that the ex- 
planation might be found in some sort of interaction 
with local structure which is not already included in the 
foreground analysis. If this is the case, then the fact 
tha t the AOE points in the direc tion of the Virgo clus- 
ter l)de Oliveira-Costa et alJ[2004fl is certainly intrig uing, 
and has led some to suggest llAbramo fc Sodrel2003|i that 
the SZ ijSunvaev fc Zeldovichlll972j) imprint of the local 
supercluster might help explain at least some of the AOE. 
Although this idea is initially attractive, the SZ is prob- 
ably 3 or 4 orders of magnitude too small to do the trick 
Ipolag et alJl200l . 

Local effects remain an attractive option, especially 
since the mass distribution is so highly anisotropic. The 
observed motions of galaxies in the local universe have 
revealed a b ulk flow in the direction o f the "Great Attrac- 
tor" region l)Lvnden-Bell et aLl fl988 K The Local Super- 
cluster, in whic h galaxies lie preferen tially in the super- 
galactic plane ijde Vaucouleursl ll9581. dominates as we 
move closer to the Earth, and of course the Milky Way 
itself dominates on still smaller scales. 

These all contribute to the velocity of the Earth rel- 
ative to the background, which gives rise to the CMB 
dipole term. The dipole is more than two orders of mag- 
nitude larger than the quadrupole and octopole, and is 
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by far the largest of the temperature anisotropics; how- 
ever, since the dipole arises from the peculiar motion of 
the observer with respect to the background, it is clearly 
of non-cosmological origin, and is therefore measured and 
subtracted from maps of the CMB. Here, we show that 
this subtraction is imperfect. Weak gravitational lensing 
by local large scale structures will coherently deform the 
initially perfect dipole, causing a leakage of power at the 
sub-percent level into other \ow-£ moments. 

We suggest, to our knowledge for the first time, that 
it is this lensing induced mixing of power from the CMB 
dipole that is the root cause of the AOE. We will in- 
troduce the basic idea, along with some relevant back- 
ground in Section |21 where we will show that the expla- 
nation is natural and qualitatively correct. In Section 
El we will present the results of a toy model simulation 
which demonstrate the effect to be of the right order of 
magnitude, and we will conclude in Section 0] with some 
discussion of the cosmological implications of the reinter- 
preted large scale CMB sky in the event that the lensed 
dipole ultimately proves to be the right answer. 

2. WEAK LENSING OF THE CMB DIPOLE 

Weak lensing of the CMB has l ong b e en a 
topic of interest to cosmologists ('e.g. iSel iakl Il996t 
Zaldarriaea fc Seliaklll99l lHull2000t IChallinor fc Lewis! 
2005j) . The effect is both simple and inescapable; all light 



which reaches us from the surface of last scattering (or 
any other source, for that matter) is deflected from its 
original path by the weak gravitational lensing interac- 
tion witli_the_jnattCT_ohstn the line of sight 
(see [Bartclmann & Schnciderl l20011 for a comprehensive 
review), and no exceptions are made for photons from 
the CMB dipole. Although the dipole owes its existence 
to the motion of the observer with respect to the back- 
ground, this makes no difference from the perspective of 
someone in the same frame of reference as the observer; 
one side of the universe is simply hotter than the other, 
and this anisotropy will be lensed. For the case of an ob- 
server comoving with the Earth, t he observed dipole term 
will be that measured by WMAP ijBennett et alJl20 03b). 
which is more than two orders of magnitude larger than 
the quadrupole term (and is by far the largest anisotropy 
in the CMB), so that even sub percent level scatter will 
strongly effect the low-£ moments. Also, because the 
dipole is coherent over the whole sky, it will couple best 
to lensing effects that are also coherent over much of the 
sky, so that local structures will be the dominant lenses. 

This cannot be accounted for by simply subtracting the 
measured dipole; lensing will scatter the initially pristine 
dipole into something that is only almost a perfect dipole, 
so that if we fit a dipole to the measured sky and subtract 
it, we are going to be stuck with a residual. It is this 
residual which we believe is a likely culprit to explain 
the AOE, as we will now explain. In the following, we 
will adopt a frame of reference that moves with the lens 
and observer, work in comoving coordinates \, set c = I, 
and assume a flat universe. 

2.1. Lensing and the CMB 

The notion that gravitational lensing by large scale 
structure will scatter pow er between l- modes of the CMB 
is certainly not a new one (|Seliak 1996), and while we will 
consider effects on much larger than usual angular scales, 



we will begin our discussion in the usual way. While en 
route to us, the photons of the primary CMB are gravi- 
tationally lensed (see e.g. IChallinor fc Lewisll2005l for a 
recent discussion of CMB lensing) by the mass distribu- 
tion along the line of sight all the way from the epoch of 
last scattering to the present, so that the photons we see 
have been displaced from their original position on the 
sky by an angle 
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where h = (9, </>) is the angular position on the sky, Vj_$ 
is the spatial gradient of the gravitational potential per- 
pendicular to the path dx of a given light ray, and x' is 
the distance to the surface of last scattering. The scat- 
tering can be well described as the angular gradient of 
an effective projected potential 4", so that if the source 
is much farther away than the lens, as is the case for us, 
the deflection is 
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where V is with respect to angular coordinates so that 
V = xVi, and the quantity inside the brackets defines 
. The observed temperature T is a simple re- mapping 
of the primordial value T 1 

T(h) = T'(n - a) (3) 

If we expand equation|3|to linear order, then the observed 
temperature is simply approximated as 

T(h) = T'(h) - VT • V* (4) 

where for our purposes VT is dipole as seen by a given 
lens. Equations |21 and 01 highlight some important points 
about CMB lensing. First, the lensing deflection falls 
off slowly, as a ~ 1/8, so that the deflection effect of 
large concentrations of mass extends much further on the 
sky than lensing induced shea r and magnification (e.g. 
iBartelmann fc Schneiderll2001|h Second, the magnitude 
of the deflection from a given overdensity is weighted as 
a ~ M/x, so that the relatively local universe will con- 
tribute substantially to the overall deflection field; points 
one and two are crucial given our need for deflections 
which are coherent over huge fractions of the sky. Third, 
since VT and V\P are proportional to the velocity and 
mass, respectively, the scatter will be maximized around 
large momentum overdensities. Finally, power which is 
scattered from the CMB dipole will be locally perpen- 
dicular to its gradient, establishing a preferred axis. 

2.2. Lensing and the integrated Sachs-Wolfe effect 

We have so far considered the gravitational interaction 
of the CMB with lenses that are comoving with the ob- 
server, and we will now generalize this result. We will 
start by considering the effect in the rest frame of the 
CMB, where it is known as the integrated Sachs- Wolfe 
(ISW) effect. We note that a number of people have 
commented to the author that the ISW is in some ways 
a more intuitive description. While we do not share this 
view, it's certainly worth pointing out the relation, dic- 
tated by the equivalence principal, between the two (in 
our opinion) equally valid descriptions. 
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Fig. 1. — The CMB te mperature quadrupole (top pa nels) and octopole (bottom panels). From left to right: the observed moments from 
WMAP (as reported by do Olivcira-Costa ct al. 2004), the moments from the simplest toy model with one attractor, and the moments 
from the 6 attractor model. Note that the color scale of the plots has been adjusted for display, spanning ±35mK for the WMAP octopole, 
±10/xK for the one attractor octopole, and ±15^tK for the rest. 



Let us consider a photon with a 4-momentum — 
Eo(l,n) traversing a potential which is changing in time, 
so that the energy a photon gains when falling into the 
well is not equal in magnitude to that lost when climbing 
out. The observed potential will be described by a met- 
ric with Q(x — (Vicns — Vobs)t), so that (working to first 
order) the observed change in energy is 

dEisw = 2d\E (Mens - Fobs) • VL$l cns (5) 

where V h s is the observer's velocity and Vj_ implies the 
direction perpendicular to the path of the photon. To 
compute the dipole lensing term, we take the dot product 
of the change in 4-momentum of the photon with the 4- 
velocity the observer, so that dE\ cns = -EoVobs • dn, and 

dE lcns = 2 d X E f bs • Vl$i c „ s (6) 

Adding equations [S] and integrating, and noting that 
E V\ cns is the dipole magnitude seen by the lens, we ar- 
rive at 

A£ = 2VT dipolo -d (7) 

3. RESULTS FROM A TOY MODEL 

We compute the lensing distortion of the CMB dipole 
in two steps. First, for a given mass distribution, we com- 
pute the deflection field a(h) using the 3-d ray-tracing 
method described in Vale fc W hite f 20031) over the whole 
sky. This method first computes the 3-d gradient of the 
gravitational potential on a 3-d grid, maps the result 
onto the points of a 3-d spherical grid, and selects the 
components which are locally perpendicular to the line 
of sight. This gives us Vj_<I>(x, n) from equation ^ so 
that a is given by a sum over x at each position on the 
spherical surface n = (6,(j>). We then generate the tem- 
perature map for the lensed dipole on the 2-d spherical 
grid by following the prescription of equation^ in which 
we compute the dipole temperature at position (h — a) 
and apply it to position h. The spherical harmonic co- 
efficients aim's of the new field are then computed, and 
the I > 1 terms are observed to no longer be zero. 



However, before we can employ this technique, we must 
first decide on a mass distribution from which to derive 
the lensing deflection. This is a challenging proposition 
given the current level of uncertainty in the local distri- 
bution of dark matter. The closest superclusters such as 
the Local Supercluster, the Centaurus Wall, the Perseus- 
Pisces chain, and the Great Attractor, all intersect the 
Milky Way "Zone of Avoidance" , so that Galactic extinc- 
tion makes a full probe of their propertie s currently prob- 
lematic (see e.g. lKraan-Kortewedl2005l for a recent dis- 
cussion). Furthermore, even if this obstacle were some- 
how overcome, our ability to probe the dark matter dis- 
tribution, using optical tracers or any other method, is 
highly uncertain. We are at present unable to confidently 
state even the source of the mass dipole moment respon- 
sible for the bulk flow of galaxies in the local universe, 
with opinions divided between the so called Great Attrac- 
tor region (centered on the ultra-massive Norma cluster 
ACO 3627), which is roughly 45/i _1 Mpc from us, and 
the Shapley Supercluster, which is 100/i _1 Mpc further. 

Because we are interested in a proof of concept anal- 
ysis, we avoid the painstaking complexities involved in 
modeling the local universe, which cannot in any case be 
done with any more certainty than current observations 
allow. Instead, we employ a toy model motivated by the 
flow of mass in the local universe. We begin by placing 
an overdensity near the Shapley concentration at a dis- 
tance of 200/i -1 Mpc, with a fwhm of 100/i _1 Mpc and a 
peak density five times the mean density of the universe, 
which is consi stent with the observed galaxy overdensity 
in the region fDrinkwa ter et aLll2004l) . The overdensity 
is sufficient to induce a peculiar velocity at the Milky 
Way of roughly 600 km/s, which increases to over 2000 
km/s near the peak of the mass; re markably, this is also 
consi stent with recent observations l)Kocevski fc Ebelind 
2005). The magnitude of the induced quadrupole and 
octopole are similar to those found in the WMAP data, 
about equal for the former and one third for the latter, 
while the orientation of these (Figure ^) is also remark- 
ably close given the simple model employed. 
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As a demonstration (Figure 0, the author has at- 
tempted to recreate the observed I — 2,3 modes "by 
hand" by placing several ovcrdcnsitics and voids at a 
distance of 200/i _1 Mpc, smoothed now on a 150/i _1 Mpc 
scale, with no net mass added. The resulting density field 
peaks at three times the mean density in the direction 
of the Shapley concentration, where the velocity of the 
flow also peaks at just over 2000 km/s, and the induced 
peculiar velocity at the Galaxy is 700 km/s. Although 
this toy is clearly not meant to be taken seriously, it is 
instructive that no particularly unrealistic mass concen- 
trations are required to reproduced the observed modes. 
It is also instructive to note that most of the signal comes 
from structures located far away from any overdensities; 
the velocity field effects a huge volume of space, which 
is typically at the mean density, and it is this extended 
region that contributes the lion's share of the effect. 

We note that it is also relatively simple to create dis- 
tributions which contribute power to alternating modes; 
for example, a dumbell configuration, with the Galaxy at 
the center, contributes power to the I = 2, 4, 6 . . . modes, 
and none at all to the odd modes. We also note that 
because power is preferentially added in the direction of 
the flow, this may fully or partially explain the north- 
south power asymmetry (Eriksc net alJl2004|) . However, 
these effects depend strongly on the specific distribution 
of mass, so we have not attempted to quantify them here. 

4. DISCUSSION 

We have introduced a novel explanation for the align- 
ment of the axes defined by the planes of the CMB 
quadrupole and octopole, with each other and with the 
CMB dipole, which has become known as the "Axis of 
Evil" . We find that the CMB dipole induced by the flow 
of local structures is scattered by gravitational interac- 
tion into the higher moments, and that although the ef- 
fect on the dipole is at the sub-percent level, it is of order 
one for the higher moments, whose measured values are 
roughly two orders of magnitude smaller than the typical 
dipole seen by local structure. We have shown that dipole 
lcnsing is a natural mechanism to align axes and create 
an asymmetric distribution of power, and have explicitly 
demonstrated that it is not out of bounds to attribute 



the entire measured CMB quadrupole to dipole lensing. 
Consequently, it is highly unlikely that dipole lensing can 
be safely ignored in attempts to measure the large angle 
power in the CMB. 

While we may have mitigated the alignment problem, 
the mechanism will on average add power to the low--£ 
moments, so that the low quadrupole anomaly is stronger 
than ever. Low-^ modes are profoundly affected by the 
influence of Dark Energy (DE) through the integrated 
Sachs- Wolfe effect. At late times, most models predict an 
epoch of equality between DE and Dark Matter (DM), 
after which DE causes an accelerated expansion. This 
expansion causes potential wells on scales not yet viri- 
alized to decay, causing a differential in the blue and 
red-shift photons experience as they traverse in and out 
of the wells. Overestimating the low-^ power will lead to 
an ovcrcstimation of this differential, and lead to erro- 
neous conclusions about the epoch of DE-DM equality, 
and therefor about dark energy parameters. 

Overestimation of the low-£ modes of the temperature 
may cause still other problems. One of the many im- 
portant observations in the WMAP data is that there 
is more than the expected amount of power in the low-£ 
modes of the cross correlation, Cj E , between E-polarized 
modes and temperature, which has been used to infer a 
higher than expected re-ionization depth l)Bennett et alJ 
2003b). If much of the power in the observed data is 
due to the dipole and lensing, then it may be that the 
measured Cj E is biased high, which would alter the pre- 
dicted ionization depth. 

Whatever their cause ultimately proves to be, the low- 
£ anomalies in the CMB temperature remain one of the 
fertile areas of research beyond the standard cosmologi- 
cal model, and just might prove to be the right window 
to look through to glimpse the physics driving the ex- 
pansion of the universe. 
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NASA. 
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